Vectorized numerical algorithms for the solution of continuum mechanics problems
Annotation
The aim of the work is to study the possibilities provided by new information technologies, object-oriented programming tools and modern operating systems for solving boundary value problems of continuum mechanics described by partial differential equations. To discretize the basic equations, we applied the method of finite differences and finite volumes, which are widely used to solve problems in the mechanics of liquids and gases. The paper considers the implementation of the finite difference methods and the finite volume method with vectorized grid structures, including access to the inner and boundary cells of the grid, as well as the features of the implementation of algorithms at singular points of the computational domain. To solve boundary value problems described by partial differential equations, we developed an approach to the construction of vectorized algorithms and considered the features of their software implementation in the MATLAB package. Vectorization in such tasks, excluding nested loops, is ensured by appropriate data organization and the use of vectorized operations. On the one hand, the developed algorithms widely use MATLAB functions designed for processing vectors and sparse matrices, and on the other hand, they are distinguished by high efficiency and computation speed, comparable to those of programs written in C/C++. The main results imply the numerical solution of a number of problems in continuum mechanics associated with the calculation of stresses in a separate body and the calculation of the field of velocity and temperature in the flow of a viscous incompressible fluid. The features of discretization of the basic equations and the implementation of the corresponding finite-difference and finite-volume algorithms are shown. The use of the MATLAB system opens up new possibilities for the formalization and implementation of finite-difference and finite-volume methods for the numerical solution of boundary value problems in continuum mechanics. Despite the fact that the capabilities of the developed algorithms are illustrated by the example of solving fairly simple problems, they admit a relatively simple generalization to more complex problems, for example, solving the Euler equations and Navier–Stokes equations. As part of the work, computational modules were prepared with user programming tools that expand the capabilities of the package and are focused on solving problems in continuum mechanics.
Keywords
Постоянный URL
Articles in current issue
- A study of a silicone film deposited on quartz glass under laser radiation
- Optical composites based on organic polymers and semiconductor pigments
- A new algorithm for the identification of sinusoidal signal frequency with constant parameters
- A study of silicon p-n structures with mono and multifacial photosensitive surfaces
- Detection of yawning in driver behavior based a convolutional neural network
- A Game Theory approach for communication security and safety assurance in cyber-physical systems with Reputation and Trust-based mechanisms
- A study of the influence of human factors on the speed of urban rail transport
- An algorithm for detecting RFID-duplicates
- Reduction of LSB detectors set with definite reliability
- Classification of objects in images with distortions based on a two-stage topological analysis
- Dimensionality reduction of the attributes using fuzzy optimized independent component analysis for a Big Data Intrusion Detection System
- An optimal swift key generation and distribution for QKD
- A study of vectorization methods for unstructured text documents in natural language according to their influence on the quality of work of various classifiers
- Recognition the emotional state based on a convolutional neural network
- Intellectualization of personnel development management in high-tech service-oriented companies
- A study of the efficiency of the magnetic compass correction system
- A new analytical model of drain current and small signal parameters for AlGaN-GaN high-electron-mobility transistors
- Imputation and system modeling of acid-base state parameters for different groups of patients
- Construction of movement trajectories for objects based on the Dubins car problem, taking into account constant external influences
- A mathematical model of an epidemic with an arbitrary law of recovery
- Simulation of the pulsed outflow of air and fine powder mixture, partially filling the discharge channel
- A comparative analysis of computational intelligence algorithms for estimation of LTE channels
- Implementation of a clinical decision support system to improve the medical data quality for hypertensive patients